Signatures of critical full counting statistics in a quantum-dot chain 
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We consider current shot noise and tiie full counting statistics in a chain of quantum dots which 
exhibits a continuous non-equilibrium phase transition as a function of the tunnel couplings of the 
chain with the electrodes. Using a combination of analytical and numerical methods, we establish 
that the full counting statistics is conventional away from the phase transition, but becomes, in a 
well-defined sense, essentially non-Gaussian on the critical line, where the current fluctuations are 
controlled by the dynamic critical exponent z. We find that signatures of the critical full counting 
statistics persist in quantum-dot chains of finite length. 
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I. INTRODUCTION 

It is by now well established that non-equilibrium cur- 
rent fluctuations in nanoscopic conductors yield much in- 
formation which is not contained in the current-voltage 
characteristics. Most prominently, the second moment of 
the non-equilibrium current fluctuations, known as shot 
noisCji is sensitive to effects of quantum statistics on the 
current flow and provides access to the charge of excita- 
tions. Experiments accessing the higher moments of the 
current fluctuations are now becoming available;^ in some 
cases including measurements of the entire full counting 
statistics.—!^ 

The full counting statistics^i^ (FCS) generalizes the 
concept of photon counting statistics in quantum op- 
tics to nanoscopic conductors and characterizes the cur- 
rent fluctuations by means of the entire distribution 
function Pt{Q) of the charge Q passing through the 
conductor during a time interval of length T. The 
FCS has been investigated in a wide variety of sys- 
tems, including but not limited to superconducting^ 
and normal-superconductor hybrid systems^ tunnel 
junctions»2, chaotic cavities^ spin-correlated systems^^ 
quantum dots in the Coulomb blockadoi^ and Kondo 
regimes;^ single-molecule junctions r^ii^ and nanoclcc- 
tromechanical systems 

A common feature of the FCS in all of these systems is 
that any cumulant ((Q")) of the FCS is proportional to 
T. This implies that in a well-defined sense, the FCS is 
essentially a Gaussian distribution with only small devia- 
tions. In fact, if we express the FCS in terms of the vari- 
able q^{Q~ {Q})/Vt which measures the fluctuations 
of the transferred charge in units of its typical magnitude 
as determined by the variance - V^, then all 

cumulants ((9")) of the rescaled FCS Pt{<i) other than 
the variance tend to zero as T — > 00. 

It is the purpose of the present paper to show that 
the FCS of a chain of quantum dots can behave in a 
fundamentally different manner, having cumulants (n > 
2) which scale as 
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FIG. 1: (Color online) (a) Schematic depiction of the 
quantum-dot chain, (b) Mean-field phase diagram of the dot 
occupations, (c) T-dependence of variance {(Q^)) both on 
(a = /3 = 0.5) and away from (a = 0.3, (i = 0.7) the critical 
line. 



in an appropriate "thermodynamic limit." This means 
that all cumulants of the FCS, when rescaled by the vari- 
ance as described above, are of the same order, implying 
that the FCS is essentially non-Gaussian. This result 
is a consequence of a continuous non- equilibrium phase 
transition occurring in the system. Its critical indices are 
known exactly and the exponent appearing in Eq. ([1]) can 
be identified with n/2z where 2: = 3/2 is the dynamic 
critical exponent. We explore the FCS of the system and 
its consequences for current noise by a combination of 
analytical and numerical techniques, paying particular 
attention to an analysis of finite-size corrections which 
would be relevant in experimental realizations. 



II. MODEL 

Consider a series of N weakly coupled metallic quan- 
tum dots as shown in Fig. [U under the following condi- 
tions: (i) The voltage between neighboring dots is suffi- 
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cicntly large and energy relaxation within the dots fast 
enough that transport is unidirectional, i.e. electrons al- 
ways tunnel in the direction of the voltage bias (say to 
the right), (ii) Due to the Coulomb blockade, the quan- 
tum dots can only switch between two neighboring charge 
states which we denote by n = and n = I. (iii) The 
tunneling rates between dots are equal to Fg, while the 
rate for entering the first dot from the left reservoir (leav- 
ing the last dot to the right reservoir) is tunable (e.g. by 
a gate electrode) and equal to a {(3). A more detailed 
discussion of these conditions including the electrostatics 
of the quantum-dot chain is relegated to the Appendix. 
When ignoring correlations between the occupations rii 
of different quantum dots within a mean-field description, 
the rate of tunneling between quantum dots i and i -I- 1 is 
given by TQni{l~ rii+i). Thus, we can describe transport 
through this quantum-dot chain by the rate equations 



To [ni-i{l - rii) - n,,(l - n^+i)] + SJi - SJ^+i (2) 



for 2 < i < {N — 1). The occupations of the first 
and last dot are similarly determined by the equations 
drii/dt = a(l — ni) — roni(l — 712) + SJi — SJ2 and 
dnpf/dt = ronjv-i(l — tin) — Pn^ + 6Jn — SJn+i- Since 
we will be interested in computing current fluctuations, 
we have already included Langevin sources SJi into the 
rate equations which account for the stochastic nature of 
the tunneling processes. In analogy with the Boltzmann- 
Langevin approach to fluctuations of the distribution 
function ; ^^1^^ the Poisson nature of the tunneling pro- 
cesses gives {6Ji{t)6Jj{t')) = roni_i(l — ni)Sij6{t — t') 
for the correlation function of the Langevin sources. 

The stationary states of Eq. ^ are well known from 
studie a^^'^" of the totally asymmetric simple exclusion 
process (TASEP) of which the quantum-dot chain under 
consideration is a particular realization. Remarkably, the 
stationary states exhibit different phases as function of 
a and (5. This is summarized in the phase diagram Fig. 
[1] which exhibits low- and high-density phases I and II, 
which are separated from a maximal current phase III 
by a continuous non-equilibrium phase transition. (The 
transition between the high- and low-density phases is 
first order.) These results are well established by exact 
solutions and simulations r^i^ but can also be obtained 
in essence at mean-field level. 

The mean-field occupation profile Ui can be derived 
from a recursion relation based on current conservation 
c ~ ni(l — m^i)^ (Here and in the following, we use 
units of time such that Fo = 1.) One finds that the oc- 
cupations Ui = a {rii = 1^/5) in the low-density (high- 
density) phase are constant in the bulk. Near the bound- 
ary of the chain, they relax to the bulk value (denoted 



by n in the following) within a distance ^ = 1/2(1 — 2n). 
The divergence of the correlation length ^ at n = 1/2, 
i.e., along the phase transition line between the high- or 
low-density phase and the maximal-current phase, signals 
the occurrence of a continuous phase transition. In anal- 
ogy with second-order phase transitions in thermal equi- 
librium, ^ can also be extracted from a linearized mean- 
field description of the density-density fluctuations of the 
model based on Eq. From these results, one also 

concludes that the respective average currents through 
the chain are equal to c = a[l — a) (c = /3(1 — (3)). 

The current fluctuations, including the entire distri- 
bution function, have been discussed for several related 
models in a number of publications i^^'^'^i^^'^^i^^'^'' These 
works focused on exact solutions for infinite chains and 
ring geometries with particular choices of initial condi- 
tions. Here, we employ an alternative field theory de- 
scription starting from (the approximate) Eq. ^ which 
builds on existing approaches to the FCS of nanoscopic 
conductors'^ and which also provides a convenient start- 
ing point to discuss finite-size systems. 



III. FIELD-THEORY DESCRIPTION OF THE 
FCS IN THE THERMODYNAMIC LIMIT 

In terms of the charge Qi passing the bond between 
quantum dots i — 1 and i, 



di[n,_i(l-n,)+(5J,], 



(3) 



the FCS (generalized to the joint distribution function of 
all the Qi) is defined by 



Prm.}) 




dt[n,_i(l-n,)-fW.] ) (4) 



In performing the average over the (5 J, we need to remem- 
ber that the occupations rii must satisfy the Boltzmann- 
Langevin equation ([2]) and are thus themselves dependent 
on the 5J . We can formally average independently over 
occupations and Langevin sources after enforcing the 
Boltzmann-Langevin equation by means of i5-functionals 
in the average. Doing so, passing to the Fourier transform 
PT{{X^}) = /n,dQ.exp(iE, XjQj)^T({Q,}), and ex- 
ponentiating the functional (5-functions, we obtain^^ 



Pt{{x^])^^ jX{[Ao.m (c'^)„ 



SJ 



(5) 



with the action 



dnj_ 
dt 



nj-i(l - rij) + nj{l - n^+i) - 5Jj + 5Jj^ 



(6) 
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The prefactor Z ensures the normahzation condition PriXxi = 0}) = 1- Moreover, we have made the counting fields 
Xj time dependent for simphcity of notation. PT({Xj}) follows by setting Xi(^) = Xi for < t < T and zero otherwise. 
Performing the Gaussian average over J J, passing to the continuum limit (with the distance between quantum dots 
set to a = 1), and integrating over the constraint field a^, wc obtain after lengthy, but straight-forward manipulations 

AxAt I xd-^dtn + ^ (dtn + - lain") [dtn + - \ (7) 



where c = n(\ — n) denotes the average current in terms 
of the average occupation n. For simplicity, we removed 
a trivial term from this action such that its Fourier trans- 
form gives a shifted FCS defined as the distribution func- 
tion of Q — cT . a quantity with zero average. In the fol- 
lowing, this shifted charge variable will be denoted by 

It is important to note that the action Eq. ([7]) is non- 
Gaussian in the occupation field n due to the presence of 
the terms involving j = n(l — n). Writing n = n -f- An, 
we have d^j = (1 — 2n)dx^n — One finds 

by power counting that the non-linear convective terms 
arc irrelevant away from the critical line well inside the 
high- and low-density phases where fi ^ 1/2, but become 
relevant at the continuous phase transition where n = 
1/2. We now turn to an analysis of the FCS based on 
the action Eq. 

Employing the action ([7]) in the absence of the non- 
linear convective terms yields a Gaussian FCS which is 
fully described by the variance of Q. Using {{Q"^)) — 
- afr In -Pt(x) I x=o, we obtain 

((O^)) =c [ ' dm' 

J-T/2 



duj 1 -k ^ 



(8) 



in terms of the drift velocity w = 1 — 2n. Performing the 
time and frequency integrations gives 



dq 1 - e-'?'^/^ cos{vqT) 



27r 



q- 



(9) 



For nonzero v, i.e., well inside the high- and low-density 
phase, the exponential factor in the numerator can be 
replaced by unity and wc obtain ((Q^)) = cvT. 

The result for ((Q^)) is strikingly different on the crit- 
ical line where v = 0. Here, one finds from Eq. Q 
(i.e., within hnearized mean-field theory) that ((Q^)) = 
C\J2T j-K. The crucial observation is that the variance no 
longer scales proportional to T. Eq. ([9]) suggests that we 
can identify the exponent of T in this relation as I/2 in 
terms of the dynamic critical exponent z which takes the 
value z = 2 within linearized mean-field theory. 

Indeed, this identification allows us to extend our re- 
sults beyond linearized mean-field theory. For w = 
the full action Eq. ([7]) is closely related to the one- 
dimensional Burgers equation 
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FIG. 2: (Color online) Full counting statistics of the totally 
asymmetric exclusion process vs a scaled charge variable Q on 
the critical line (main panel) [a = /3 = 0.5], where the FCS 
is essentially non-Gaussian and far from the phase transition 
(inset) [of = 0.3, /? = 0.7], where the scaled FCS converges to 
a Gaussian for large T . 



for a field (\)[x, t) driven by a Langevin source ri[x, t) with 
correlator {r]{x, t)rj{x' ,t')) ~ 5{x — x')5{t — t'). This sug- 
gests that the true variance of the quantum-dot chain 
obeys ((Q^)) ~ T^/^, where z is the exact dynamic crit- 
ical exponent z = 3/2 of the Burgers equation.— This is 
confirmed by our numerical simulations as shown in Fig. 
[TJ:. 

Within linearized mean-field theory, the FCS satisfies 
the scaling relation Pt{Q) = T-^/^f{Q/T'^/^) in terms of 
a scaling function f{x). This can be seen directly from 
the linearized version of the action Eq. ([7]) with w = 
by noting that it remains invariant under the rescalings 
t t/T, X -> x/ri/2, n uT^^'^, and x ^ xT^^'^- In 
view of our results for the variance, one may then expect 
that the exact FCS satisfies the scaling relation 
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(10) 



with z ~ 3/2. This is indeed nicely confirmed by our nu- 
merical simulations as shown in Fig. [21 This scaling rela- 
tion is at the heart of the anomalous FCS of the quantum- 
dot chain under considerations. Indeed, the result for the 
cumulants in Eq. ([T|) is an immediate consequence of this 
relation. 
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FIG. 3: (Color online) Current noise as function of fre- 
quency CO (a) well away from the critical line (with parameters 
a = 0.2, j3 = 1 — a) and (b) at (near) the critical line with pa- 
rameters a = 0.5 (q — 0.425) and /3 = 1 — a. The simulation 
data are compared with the predictions of linearized mean- 



field theory (Eq. [T3)l and the critical scaling S{uj) 



The assignment of symbols and lines is given in the legends 
within the figures. 



From the point of view of current experiments, it is 
also interesting to calculate the power spectrum S(uj) = 
2 J dte''^''{SI{t)SI{t + t)) of the current I{t). Far from 
the critical line, we readily obtain 



{uj - vqY + g4/4 



(12) 



within linearized mean-field theory. Evaluating the sum 
over q yields 

S{lo) = c[v2Rc(w2 - 2iw)-i/2 + Re(t;2 - 2ic<j)i/2] (13) 

which leads to a constant noise power S{uj) = 2cv for 
small Lu. Sufficiently far from the critical line, this result 
is in good agreement with our simulations, sec Fig.[3fa). 

On the critical line, linearized mean-field theory pre- 
dicts 



(14) 



The particular value of the exponent in this relation is an 
artifact of linearized mean field theory. The true scaling 



of S{uj) at small frequencies can be extracted from the 
relation ((Q^)) = J{duj/Truj^)S{uj)sm'^{ujT/2) combined 
with Eq. Il]) from which we read off that S{lu) cx lu^^^^^ . 
This is confirmed nicely by our simulations, as shown in 
Fig.EJb). 

In the vicinity of the critical line, i.e., for a small but 
finite drift velocity v, the noise power S{uj) crosses over 
between the linearized mean-field theory results for small 
frequencies and the critical behavior for larger frequen- 
cies, see Fig. ^h). This can be understood in terms of 
the correlation length ^ = l/2v. Once uj > v"^, excita- 
tions do not drift far enough during a time lu~^ to explore 
the finite correlation length of the system, v/oj < ^. As 
a result, linearized mean- field behavior applies only for 
u! < while the system appears critical for lo > . 



IV. FINITE-SIZE CORRECTIONS 

Any experimental realization would contain only a fi- 
nite number TV of quantum dots. In this section, we dis- 
cuss to which degree signatures of the critical full count- 
ing statistics persist in finite-size systems. For finite A^, 
one expects that the full counting statistics obeys Eq. 
([TT|) up to times T ~ . Within linearized mean- field 
theory, this can be obtained by the following argument. 
In a finite-length system, the upper limit for the corre- 
lation length ^ = l/2f is ^ = A", from which we can ex- 
tract an effective drift velocity v given by v^s = 1/{2N). 
Thus, we expect finite-size effects to become relevant once 
DoffT = A^ which immediately yields the above estimate 
with z = 2. Thus, the anomalous scaling persists up to 
values of T which are large multiples of the microscopic 
tunneling rate I/Fq, even for moderate values of TV. 

It turns out that the signatures of the critical elec- 
tron dynamics remain most pronounced in the variance 
and the noise power. This is shown in Fig. lUJa) which 
exhibits the variance ((Q^)) as function of T for small 
system sizes. It is apparent from the plots that a sub- 
linear T dependence is observed for small T, for chains 
as short as TV = 4. The dependence is close to the non- 
linear behavior expected in the infinite system when T is 
not too large, and crosses over to linear behavior only for 
larger T in accordance with the discussion in the previ- 
ous pargraph. Indeed, we can extract the maximal T up 
to which critical behavior persists from ((Q^)) by plot- 
ting ((Q2))/r2/3 vs T, see the inset of Fig. g^b). In 
this plot, critical full counting statistics is indicated by 
a constant while deviations from the constant at large T 
indicate significant finite-size corrections. Extracting a 
Tmax for different TV and plotting Tmax vs TV nicely fits 
the expected T^ax ~ TV^^^ dependence with the exact 
exponent z = 3/2, see Fig.|Ub). 

Experimentally, it may be more accessible to study 
the noise power as a function of frequency. At small 
frequencies uj < T^~^^ ^ ^/N", the behavior should be 
noncritical while the critical scaling S{uj) ^ uj^^^ would 
persist at larger frequencies. This expectation is borne 
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FIG. 4: (Color online) (a) Dependence of variance {(Q^)) on 
time T for small system sizes A*' = 4 and 10. For comparison, 
a system with A'^ = 1000 is also shown, (b) Crossover time 
Tmax up to which critical scaling persists as function of system 
size iV. The points are extracted by plotting {{Q^))/T^^^ vs T 
at each A'' and identifying the point where the curve deviates 
from a constant, as illustrated in the inset for a particular 
value of A'^. (Note that the linear fit of the curves for T > 
Tmax merely serves to identify Tmax and does not reflect a 
theoretically expected result.) 



out by our numerical results, see Fig.|5l 

One also expects the transition to smear over a certain 
range of tunneling rates a and (3 when N is finite, im- 
plying that in finite-size systems signatures of the critical 
electron dynamics should be visible even away from the 
nominal transition line. The region over which the transi- 
tion becomes smeared can be estimated as follows. As ar- 
gued above, critical behavior is visible when wTmax ^ N ^ 
i.e., for V < N^~^. Using that v — 1 — 2n, we find a criti- 
cal region \n — 1/2| < 1/{2N^/^). This estimate is consis- 
tent with our numerical results as shown in Fig. jB] This 
figure shows that the dependence of the variance {{Q^)) 
on T exhibits sublinear behavior in the entire critical re- 
gion. Moreover, it is demonstrated in the inset that the 
dependences are rough power laws with exponents close 
to 2 /3 throughout the critical region, while the exponent 
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FIG. 5: (Color online) Noise power as function of frequency 
for quantum-dot chains with A'^ = 10, 20 and 50. Note that 
one recovers the critical scaling for sufficiently large frequency 
for all system sizes. 
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FIG. 6: (Color online) Variance {{Q^)) as function of T for 
finite A' = 10, at various distances from the nominal criti- 
cal line. Sublinear behavior is observed in the entire critical 
region, as explained in the text. Within the critical region, 
the dependence be approximately fitted to a power law, with 
exponents close to 2/3, as illustrated in the inset. 



approaches unity outside of the critical region. 



V. CONCLUSIONS 

We have shown that a chain of quantum dots may re- 
alize a totally asymmetric exclusion process and thus ex- 
hibit a non-equilibrium phase transition as function of 
the tunnel couplings to the leads which results in critical 
full counting statistics of the current fluctuations. We 
find that the non-equilibrium phase transition controls 
the dynamics over a wide range of time scales, even in 
quantum-dot chains of finite length. Our minimal model 
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FIG. 7: Equivalent circuit of the quantum-dot chain. White 
(black) dots indicate quantum dots (electrodes). 

assumes that the hopping rates between quantum dots 
are all equal. Due to the exponential sensitivity of tun- 
neling, this would presumably be exceedingly difficult to 
realize in a lateral arrangement of quantum dots. In- 
stead, it appears more promising to employ a vertical 
setup)^ in which the tunneling barriers and the quantum 
dots can be formed by a well-defined number of mono- 
layers. We expect that including weak backscattering or 
rare double occupation of a dot would leave our results 
qualitatively unchanged. Nevertheless, it would be inter- 
esting to study their consequences as well as the influence 
of variations in the tunneling rates in more detail. 

It would also be interesting to identify other systems 
which exhibit a critical full counting statistics, especially 
systems whose dynamics is more directly controlled by 
quantum mechanics. An intriguing possibility is trans- 
port near a quantum phase transition. 
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APPENDIX A: ELECTROSTATICS OF A 
QUANTUM-DOT CHAIN 

In this appendix, we briefly discuss the electrostatics 
of a quantum-dot chain. We consider metallic quantum 



dots with a continuous spectrum so that the tunnel rates 
can be assumed to be a linear function of the energy 
gained in a tunneling process. To derive this energy gain, 
we model the quantum-dot chain by the equivalent circuit 
shown in Fig. [7] which is described by the capacitance 
matrix C defined as 




where the index d(e) enumerates the quantum dots (elec- 
trodes). The charges and the potentials on the dots and 
electrodes are denoted by q and v, respectively 

Up to irrelevant constants, the energy of the system 
takes the form 



U = -Vd^CddVd. (A2) 

In an infinite chain, the energy difference At/; induced 
by a tunneling event from the zth to the (i -t- l)th dot 
becomes 

AC/, - -V, + l±^h±L^ (A3) 

, C ( ni+2 - rii+i + rtj - nj_i - 1 \ 

^c; [ '^"j ' 

to first order in the weak coupling limit C ^ Cg. Here, 
{rii} denote the occupation numbers of the initial state 
and we use units with e = 1. The weak coupling limit 
C <^ Cg ensures that the contribution to the energy aris- 
ing from interdot interactions is small compared to the 
contribution of the gate capacitances. According to Eq. 
(|A3p we conclude that double occupation and backscat- 
tering is strongly suppressed as long as fc^T <C Vq <C 
e^/Cg while the tunneling rates are approximately equal 
as long as C <C Cg. 
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